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Abstract 

It is known that a linear hamiltonian system has too many invariant measures, thus the problem of 
convergence to Gibbs measure has no sense. We consider linear hamiltonian systems of arbitrary finite 
dimension and prove that, under the condition that one distinguished coordinate is subjected to dissipation 
and white noise, then, for almost any hamiltonians and almost any initial conditions, there exists the 
unique limiting distribution. Moreover, this distribution is Gibbsian with the temperature depending on the 
dissipation and of the variance of the white noise. 



1 Main Results 

Consider the phase space 

L = L2N = = {V- = (g,p) : q={qi,..., ^at), p = (pi, . . . ,Pn) & R^} 



with the scalar product 



N 



The space L is the direct sum L — l''^ ® of orthogonal coordinate space and momentum space with the 
induced scalar products {q, q')i and (p,_p')i correspondingly. We are most interested in the case of large N, but 
we do not use it in the present paper. 

We shall study the following system of 2N stochastic differential equations (fc = 1, . . . , N) 

dqk = Pk dt, 

dPk = j:fL,(^hV{k,l)qi-D{k,l)pi)dt + B{k,l)dwu), 

where V = {V{i,j)) is a positive definite (N x iV)-matrix, D — {D{i,j)) is the non-negative definite symmetric 
(N X A^)-matrix, B ~ {B{k,l)) is an arbitrary real matrix, wtj, I ~ l,...,N, are the standard brownian 
processes, independent in I. 

If _D = _B = 0, then the system is a linear hamiltonian system with the quadratic hamiltonian 

N N 

HW^T + U, T^-Y,pI U^-J2vit,j)q,q,. (2) 

i—l — 1 

We will consider systems where the matrices B and D arc as follows (for some n = 1, . . . , N) 

D = aA„, B = (tA„, 

where A„ = {5i^j5i^n) is the matrix with one diagonal element equal to 1 and other elements are all zero. Note 
that in this case the system will be subjected to only one white noise dwt,n- We assume that the index n is 
fixed, and write for shortness dwt = dwt.n- Further on, unless otherwise stated, we assume that a > 0, ti > 0. 
System (P) can be rewritten in the matrix form 

dip = Alp dt + agn dwt , (3) 
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where 

A-( ^ 
^~[-V -D 

E is the unit (N x 7V)-niatrix, 5„ = (0, e„) e zjj''', and e„ = (0, . . . , 0, 1, 0, . . . , 0) G is n-th standard basis 
vector. The solution of the latter equation with arbitrary initial vector ip{0) is uniquely defined and can be 
written as (for example, see section 12.4 of [8]) 

i^it) = e*^ (^a ^ e-'^gn dw, + ^-(0) 
Introduce the set 

L^^{il>eL: iJ(e*'^V) ^0, t oo} C L. 
We will need the following results. 

Lemma 1. The set i_ is a linear subspace of the space L, and moreover L- — {{q,p) £ L : q G ly, p G Iv} 
where ly -is the subspace of M.^ , generated by the vectors V'^en, k = 0,1, . . . , in particular gn G Moreover, 
L_ and its orthogonal complement, denoted further by Lq, are invariant with respect to the operator A. 

All assertions of this lemma have been proven in jj] . 

By this lemma any initial vector ip{0) can be uniquely decomposed as 

V'(O) = Vo + V'-, ^JoeLo, V- e L_. 

Then the solution ip{t) of the stochastic equation with initial vector '(/'(O) for any t > can be decomposed 
as 

m = ^^°Ht)+4'^~\t), (4) 

where V'*-°-'(Oi V'^^H^) satisfy the equations 

ij^"\t)^Aij^-°\t), di^'^-\t) = A^p'^-Ut + agndwt 

with the initial conditions ^''^'(O) = "00, '(/''"^(O) — ip~ correspondingly. In fact, the sum of these equations give 
the equation ([3]). 

By Lemma 1, the function ^^-^^(t) G Lq, t G [0, oo), is deterministic and can be written as 

i^(o)(i) = e^*^(°H0), (5) 
and 4'^~\t) is a gaussian random process with values in L_ (as .g„ G L_). 
Theorem 1. For any ip{0) the convergence in distribution 

takes place, where ^ G and its distribution is absolutely continuous with respect to Lebesgue measure on L_ 
( defined by the euclidean structure ), and has the following density with respect to this measure 

p^(V) = |exp(^-^i/(V)), VgL^. (6) 
The limit of the mean energy is: 

2 

lim EiJ(^/;(")(t)) = — dimL_. 
t->-+oo 4a 

Thus, the action of the random force and dissipation on one particle only garanties convergnce to the 
invariant Gibbs measure with the temperature depending on a and a. 

The first assertion of the next theorem shows that convergence to Gibbs distribution is a typical property of 
linear hamiltonian systems with dissipation and random force. The second assertion shows that the dissipative 
term is necessary for this convergence. 

For given N denote Hjv the smooth manifold of all possible hamiltonians _ff as in ([2]) , that is the smooth 
manifold of all positive definite (A^ x A^)-matrices V . Let pL be an arbitrary absolutely continuous probability 
measure on Hjv, and let H^'' be the set of all hamiltonians of H^v, for which the dimension of Lq is greater 
than zero. 



Theorem 2. 1. The measure fj. of the subset H]^ is zero. 2. If a — 0, then for any initial condition V'(O) we 
have 

2 

Eff(^(t)) = ^t + 0(1). 

Note that for more restricted (physical) classes of hamiltonians the property dim = is not typical (see 
[3] in this respect). 

In this short note we restrict ourselves to the most interesting case of one distinguished particle, which 
shows that even the minimal introducing of stochasticity to the system garanties the convergence to the physical 
equilibrium. Note however, that most results can be generalized to arbitrary matrices D and B. Similar systems, 
mainly one-dimensional (one of the goal was to justify the Fourier law of heat conduction) were considered in 
the 1960-70 in the series of papers by J. Lebowitz and colleagues (see [I], [2] and references therein). 



2 Proofs 



Proof of theorem 1. Assume first that dimio = 0. In this case the spectrum of the matrix A belongs to the 
left half-plane. 

Let us prove the convergence first. By the latter assumption, we can consider only the process ^p^~\t), 
admitting the following decomposition 

■ip(-^t) ^lP'^9\t)+lP'^'^^t), 

By definition of i- the function ijM\t) tends to zero if t — > +oo. That is why it is sufficient to prove the 
convergence of ■>j/s'>{t). Denote C{t) ~ (E {'0,p'' (i)V'j^'' (0}) the covariance matrix. Then 

C{t) = E {V'(^) m/^Yit)} = a'e'^ E I £ e'^^^n dw, £ sje"^^" dw^ 

where ^ denotes transposition. Using the Ito isometry [7], we get 

C{t) = a'e'^ f e-^^5n5je-^^" dse'^\ (7) 
Jo 

Let us calculate the integral in the last formula by finding the matrix U , not depending on time and such that 



t 

-sA T -sA"^ 1 -tArr -tA^ TT /n\ 

6 dnQn^ ds = e Ue —U. (8) 



Differentiation ([8]) in < shows that the exponents cancel and we get 

AU + UA^ = -g^gl 

This equation with respect to U has the unique solution, as the spectrum of A lies in the left half-plane, see 
section 4.4 of [5]. It is easy to check that the solution is the following matrix 

Thus, from © and dH]) we get C{t) = a'^{U - e*^f7e*^^). As the spectrum of A lies in the left half-plane, 
then limt_j.+oo C{t) — a^U. That is why the following limit in distribution 

^ lim V^^n^) = lim V-^^n*), 

exists and is the gaussian vector with zero mean and covariance matrix a^U. 

Now prove the last assertion of theorem 1. From © and from positive definiteness of V, it follows that 
U is non-degenerate. Thus the distribution of ^ has the density Ps^itp) with respect to the standard Lebesgue 
measure dqi ■ ■ ■ dpN on L. As the matrix {2aU)~^ defines the quadratic form 77, then 



Remind that H{ip) — {Qi/j,ip)i/2, where the (2A^ x 2A^)-niatrix Q — {qij) is defined by the equality of the 
vectors Q{q,p) — {Vq,p). Then for the mean energy we get 



^ 2N 

hm EHi^^-\t))=-EHiO = - ^ g,,,c,,, = — tr (QC/) - — dimi. 



Thus for the case Lq — {0} theorem 1 is proved. 

Consider now the case Lq ^ {0}. Let ui, . . . , be an arbitrary orthonormal basis of the subspace ly with 
the only restriction that vi = e„. Using it, define the orthonormal basis in L_ as follows 

hi^^^{vk,0), hf = {Q,vtd, k = l,...,d. 

The coordinates on L_ in this basis we denote V'' = {ip'i, ■ • ■ , V'm)- these coordinates our equation on L_ can 
be written as 

d^j' = A'^p' dt + g[dwt, 

where A' = _^y, ^, E is the unit matrix of the order d, V , D' and g'l are the matrices of the operators 

V ^ D and the distinguished vector, in the new coordinates, correspondingly. 

The latter equation looks like the main equation ([3]). The operator A' is the restriction of the operator A 
onto the subspace L_; let i_ = Lq L'_ be the correponding decomposition of i_ in A' . Then if, = and we 
can apply the assertions proven for the case io = 0. Note that the quadratic form H' on i_, generated by V' , 
coincides with the restriction of the quadratic form H on the subspace L_. This proves Theorem 1 completely. 

Proof of theorem 2. If a = 0, then the energy conservation law gives L_ = 0. Decompose the solution %l]{t) 

as 

where il}h{t)-\s the solution of the homogeneous equation with initial condition "0(0)1 and i^i{t) is the solution 
of the inhomogeneous equation with zero initial conditions. As H{'il){t)) has the norm properties, then 

The latter equality holds as a = 0, and thus the energy is conserved. Let us find EiJ(V'/(i)). Denote i)i{t) — 
{q''^^{t),p^^\t)). The solution with zero initial conditions is similar to the known (see for example [6^) formula 
for ordinary differential equations 

q{t) = a {Wy^ [ sm{Vv{t- s))endws, 
Jo 

p{t) = a cos{Vv{t — s))endws. 
Jo 

Using again the Ito isometry we get 

ET = ^Eipit),p{t)), = ^E{p^p} = ^ l^elcos^VV{t~s))e^ds, 
E [/ = y E iVq{t), q{t))i = ^ Eil^Vq} = y ^ ^^^^ (^(^ " ^))e« ds, 
EHiMt)) = E{T+U} = ^ j\le^ds^^t. 

We have thus proved the second assertion of theorem 2. 

Note that Hjv is a smooth manifold. Define the matrix so that its fc-th column is equal to the vector 
V^en- Then 

H^^ = {V : dim V < N} = [V : det(E(T/)) = O}. 

Note that det(S(V^)) ^ for matrices V E Hjv with simple spectrum and eigenvalue basis vi, . . . , vn, with the 
property that (vk,en)i 7^ for all A; = 1, . . . , A^, see [4l. That is why the function det{'S{V)) is not identically 
zero on the manifold Hn- Thus, H^"* is the set of zeros of the polynomial on Hn- It follows that its dimension 
is less than the dimension of all Hjv, and its measure fi equals zero. 
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